Different material models for ChamberSphere - #243
Conversation
…ged files to add this new block type to be used in a circuit model
…n the maximum number of iterations is reached
…ginal test case just with new material
…based on max./min. volume and max. pressure in pv loop from Sharifi et al. 2024
…l chamber with exponential material
…erial coice in input file. Changed ChamberSphere to accomodate different materials. Changed other files to include the new ChamberSphere and SphereMaterial structure. Changed input files to use the materials.
…ere_expmat test case)
…al parameters now. All tests pass now, but not yet included test cases for closed loop and simple expmat test case
…re is tested in closed loop test case)
![claude[bot]](https://avatars.githubusercontent.com/in/1236702?s=60&v=4){kind=small}
mrp089
left a comment
mrp089
left a comment
There was a problem hiding this comment.
Thanks, @KatrinKoesler! A view details on splitting up passive+visco and smaller readability things.
| "C1": 0.11, | ||
| "C2": 1.9e3, | ||
| "C3": 0.11, | ||
| "eta": 250.0 |
There was a problem hiding this comment.
use different sets of C0/C1 and C2/C3 for testing
| * | ||
| * Implements: | ||
| * \f[ | ||
| * f = 4(1 - C^{-3})(W_1 + C W_2) + \eta \dot{C}(1 + 2 C^{-6}) |
There was a problem hiding this comment.
Should the viscosity be separate from this? Technically, all passive materials can be combined with a viscous contribution. Since we only have one viscous model, the parameter eta could go into the sphere itself, and the viscous stress is additive to all passive materials.
| }; | ||
|
|
||
| /** | ||
| * @brief Mooney–Rivlin (neo-Hookean) material |
There was a problem hiding this comment.
The (neo-Hookean) parantheses can be misleading
| * * `steepness` - Activation steepness parameter \f$\gamma\f$ | ||
| * | ||
| * ### Usage in json configuration file | ||
| * ### Usage in json configuration file (with example material) |
There was a problem hiding this comment.
Remove all material models from ChamberSphere and use only the stress S. Reference SphereMaterial.
| const double C3 = params_.at("C3"); | ||
| const double eta = params_.at("eta"); | ||
|
|
||
| SphericalStressResult res; |
There was a problem hiding this comment.
Simplify the equations for readability? Some terms can be pre-computed
2 participants
Different material models for ChamberSphere
Current situation
As highlighted in #240 currently the ChamberSphere has a Mooney-Rivlin material model and to change to a new material model, a new ChamberSphere would have to be implemented.
Release Notes
Testing
All tests pass
Code of Conduct & Contributing Guidelines